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Abstract 

This paper concerns the asymptotic ground state properties of heavy 
atoms in strong, homogeneous magnetic fields. In the limit when the nu- 
clear charge Z tends to oo with the magnetic field B satisfying B 3> Z 4 / 3 
all the electrons are confined to the lowest Landau band. We consider 
here an energy functional, whose variable is a sequence of one-dimensional 
density matrices corresponding to different angular momentum functions 
in the lowest Landau band. We study this functional in detail and de- 
rive various interesting properties, which are compared with the density 
matrix (DM) theory introduced by Lieb, Solovej and Yngvason. In con- 
trast to the DM theory the variable perpendicular to the field is replaced 
by the discrete angular momentum quantum numbers. Hence we call the 
new functional a discrete density matrix (DDM) functional. We relate this 
DDM theory to the lowest Landau band quantum mechanics and show 
that it reproduces correctly the ground state energy apart from errors due 
to the indirect part of the Coulomb interaction energy. 



1 Introduction 

The ground state properties of atoms in strong magnetic fields have been the 
subject of intensive mathematical studies during the last decade. This paper 



is based on the comprehensive work of Lieb, Solovej and Yngvason ||LSY94a , 
|LS Y94b| ], which we refer to for an extensive list of references concerning the 
history of this subject. 
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The starting point of our investigation is the Pauli Hamiltonian for an atom 
with iV electrons and nuclear charge Z in a homogeneous magnetic field B = 
(0, 0, B) with vector potential A(x) = |B x x, 

H= E fkrHV 3 + A(x 3 ))] 2 4l+ E J^r-r ^ 



±<j<N y l<i<j<N 1 - 71 

which acts on the Hilbert space /\ 1<j<N C 2 (R 3 ; C 2 ) of antisymmetric spinor- 
valued wave functions. Here er denotes the usual Pauli spin matrices. 

The units are chosen such that h = 2m e = e = 1, so the unit of energy- 
is four times the Rydberg energy. The magnetic field B is measured in units 

2 3 

Bq = m ^3 c = 2.35 • 10 9 GauB, the field strength for which the cyclotron radius 
(-b = (hc/(eB)y/ 2 equals the Bohr radius a = h 2 /(m e e 2 ). 
The ground state energy is defined as 

E Q (N, Z, B) = inf{(tf , HV) : G domain H, *) = 1}, (1.2) 

and if there is a ground state wave function \1>, the corresponding ground state 
density is given by 

p Q (x) =p*(x) =iV E / l^(x,s 1 ,...,x 7V ,s J v)| 2 dx 2 ...dx 7V . (1.3) 

Recall that the spectrum of the free Pauli Hamiltonian on £ 2 (M 3 ; C 2 ) for one 
electron in the magnetic field B, 

H A = [<7-HV + A(x))] 2 , (1.4) 

is given by 

p 2 z + 2uB ^ = 0,1,2,..., Pz eR. (1.5) 

The projector n onto the lowest Landau band, v — 0, is represented by the 
kernel 

H (x,x') = |Uxp |i(x ± x x' ± ) ■ B - I(x ± - x'J 2 ^} *(* - ^ (1-6) 

where x^ and z are the components of x perpendicular and parallel to the mag- 
netic field, and Pi denotes the projection onto the spin-down (s = —1/2) com- 
ponent. With the decomposition C 2 (R 3 , dx; C 2 ) = C 2 (R 2 , dx ± ) <g> C 2 (R, dz) <g> C 2 
it can be written as 

n = El^)^™l® J ® p i' ( L7 ) 

m>0 

where (f> m denotes the function in the lowest Landau band with angular momen- 
tum — m < 0, i.e., using polar coordinates (r, ip), 

W 1 (Br 2 \ mt2 _..„,„ _ R , 2 , 




l ( x± ) = x l — —[-L-) e-^e- BH / 4 . (1. 
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The projector onto the subspace of Ai<j<iv £ 2 (^ 3 ; C 2 ) where all the electrons 
are in the lowest Landau band is the iV-th tensorial power of Ho and will be 
denoted by IT^. The ground state energy of electrons restricted to the lowest 
Landau band is defined as 

£ c Q onf (N, Z, B) = KHU^). (1.9) 

The assertion that for B ^> Z 4 ^ the electrons are to the leading order confined 
to the lowest Landau band is confirmed by the following theorem. 



1.1. THEOREM (Lowest Landau band confinement). f ^LSY94^ , Theo 



rem 1.2) For any fixed X = N/Z there is a 8(x) with 5(x) — > as x — > oo such 
that 

£lf >E Q > £ c Q onf (1 + 5(B/Z^)) . (1.10) 

We now define the functional we are considering. We are only interested in 
B ^> Z 4 / 3 , so we can restrict ourselves to considering all particles in the lowest 
Landau band. Recalling that 

U H A U = \<l>m){<t>m\ ® i-dl) <8> P b (1.11) 

m>0 

it is natural to define the following discrete density matrix functional (DDM) as 
= ( Tr [-^ r «J - Z I V m {z)p m {z)dz\ + D( Pl p), (1.12) 

where ^ 

D(p,p) = -^2 / Vm,n( z ~ z')p m {z)p n {z')dzdz', (1.13) 

m,n 

and the potentials V m and V m ,n are given by 

1 



V m {z) = J ^|0 m (x ± )| 2 dx ± , (1.14) 
VmAz _ zl) = f lI(xx)IV„(xl)|- w . (L15) 



Here T is a sequence of density matrices acting on £ 2 (IR, dz), 

r = (r m ) mgNo , (i-i6) 

with corresponding densities p = (p m ) m , p m (z) = T m (z,z). Note that 
depends on B via the potentials V m and V m>n . 

This functional is defined for all V with the properties: 
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(i) 

(ii) 



< T m < I for all m e N , 



J]Tr[(l-^)r m ]<oo. 

m 

The corresponding energy is given by 



1.17) 

;i.i8) 



E DDM (N, Z, B) = inf I E% U Z M [T] 



< N 



1.19) 



(We do not require that A" is an integer.) As we will show, E DDM correctly 
reproduces the confined ground state energy E® ai apart from errors due to the 
indirect part of the Coulomb interaction energy: 

1.2. THEOREM (Relation of E DDM and E%J). For some constant c A de- 
pending only on A = N/Z 



with 



> E^ on{ (N, Z, B) - E UUM (N, Z, B) > -R L 



R L = c x mm{Z 17 / 15 B 2 / 5 ,Z 8 / 3 (l + [\n(B / Z 3 )] 2 )} 



(1.20) 



1.21) 



E^ n{ is of order 



The proof will be given in Sections |5| and |6|. Note that 
mm{Z 7 / 3 (B/Z 4 / 3 ) 2 / 5 , Z 3 (l + [\n(B / Z 3 )} 2 )} , so R L is really of lower order. 
The functional flll^) 

is in fact a reduced Hartree-Fock approximation (in the 
sense of |[S91|| ) for the quantum mechanical many body problem, i.e. Hartree- 
Fock theory with the exchange term dropped. The upper bound in Theorem |1.2| 
holds quite generally for (reduced) Hartree-Fock approximations to many body 
quantum mechanics with positive two-body interactions. 

Hartree and Hartree-Fock approximations to H with restriction to the lowest 
Landau level were probably for the first time considered in [|CLR70 | and [ CJR73 1 . 
They have been studied numerically in JNTjKSBl, |JNKL87[|. 



In this paper we present a rigorous mathematical treatment of the func- 
tional (|1.12p . In Section ^| we show that there exists a unique solution to the 
minimization problem in ( |1.19| ). The corresponding minimizer is composed of 



eigenf unctions of one-dimensional effective mean-field Hamiltonians depending on 
m. Moreover, the superharmonicity of the effective mean-field potential implies 
monotonicity and concavity of the eigenvalues in m, which amounts to "filling 
the lowest angular momentum channels". This fact is important for numerical 
treatments of the model, for it means that at most the lowest angular momenta 
have to be considered. For B/Z 3 large enough, we will show that each angular 
momentum channel is occupied by at most one particle. In Section ^ we estimate 
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the maximum number of electrons that can be bound to the nucleus. We use 



Lieb's strategy to derive an upper bound analogous to ||S00|| . 

The DDM theory can also be considered as a discrete analogue of the DM 
functional introduced in [|LSY94a|| . To express this analogy we will recall its 



definition and main properties in the next section. 

2 Comparison with the DM functional 

In |[LS Y94a|j Lieb, Solovej and Yngvason defined a density matrix (DM) functional 



as 

'DM r-pl _ / rp„ _ r a 2 r 1J„ ry I I 1-1 



£ UM [T] = / Tr £2(R) [-^r x Jdx ± - Z / |x|-Vr(x) + D(p r ,p r ). (2.1) 
7r2 J 

Its variable is an operator valued function 

r :x± ^r x± , (2.2) 

where T x± is a density matrix on £ 2 (M), given by a kernel T x± (z, z') and satisfying 

< T x± < (B/2n)I (2.3) 
as an operator on £ 2 (IR). Here we denote 

p r (x) = F x± (z,z) and D(f,g) = I J ^^ dxdx'. (2.4) 

The energy 

E DM (N, Z, B) = inf {£ DM [r] | T satisfies (2.3) and J p r < N} (2.5) 

turns out to be asymptotically equal to the confined quantum mechanical ground 
state energy E® nf in the following precise sense: 



2.1. THEOREM (Relation of E DM and E^). ( ^STU^jj , Sect. 5 and 



For some constants c\ and c' x , depending only on X = N/Z , 

Ru > E cc. n {( N ' Z > B ) ~ £ DM (^ Z ' B ) > (2.6) 
with Rl given in and 



Ru = c^mmiZ^B^^Z^ll + ln(Z) + (ln(5/Z 3 )) 2 ] 5 / 6 }. (2.7) 
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The DM energy fulfills the simple scaling relation 

E DM (N,Z,B) = Z 3 E DM (X,l,r]), (2.8) 

where we introduced the parameters A = N/Z and 77 = B/Z 3 . In the limit 
rj — > 00, _E DM /(ln?7) 2 converges to the so-called hyper-strong (HS) energy E KS (\), 
which is the ground state energy of the functional 



^ HS [p] = ^(^VpR) 2 ^-p(0) + i j P {zfdz 



(2.9) 



under the condition J p(z)dz < A. The corresponding minimizer can be given 
explicitly, namely 

P [Z) ~ (4sinh[(2-A)|^|/4 + c(A)]) 2 ' 

tanhc(A) = (2 - A)/2, 
p HS {z) = 2(2+\z\)- 2 for A>2. (2.10) 

We now discuss the relation of £#£ M and S DM . The functional £j^ M is the 
restriction of £ DM to density matrices of the form 

r xx ( z ,z') = ^|0 m ( x± )| 2 r m (z,/). (2.11) 

rn 

Therefore it is clear that 

E DBM (N, Z, B) > E DM (N, Z, B). (2.12) 
In the limit N 00 they coincide, i.e. we will show that 

for all fixed A and r\. In fact, fl2.13|) follows immediately from Theorem |1.2| and 
Theorem 271 , but it could of course also be shown directly without referring to 
the relation to E® on{ . Note that in contrast to DM, the DDM energy depends non- 
trivially on all the three parameters N, Z and B, like the original QM problem. 



3 Properties of £™ 

3.1. THEOREM (Existence of a minimizer). For each N > 0, B > and 
Z > there exists a minimizer r DDM forS^ 1 ^ 1 under the condition X] m Tr[r m ] < 
N. 
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Proof. Let be a minimizing sequence for ^ M under the normalization con- 
dition ^ m Tr[r m ] < N, with corresponding densities = (pm) m - Using the 
one-dimensional Lieb-Thirring inequality |[LT76|| we have 



J |p«| 3 <const.Tr[-a 2 2 r«]. (3.1) 
Moreover (cf. ||LSY94a|| , Eq. (4.5b)), 



dJ Pm(z) 



2 



dz 



d^<Tr[-^ 2 r«]. (3.2) 



Since the potential energy is relatively bounded with respect to the kinetic energy 
(compare Thm. 2.2 in ||LSY94a|1 ), the right hand sides of ( |3.1| ) and (|3.2j ) are 



uniformly bounded. Hence the sequence pm is bounded in £ 3 (M, dz) fl£ 1 (M, dz), 
and J pm is bounded in H 1 (R, dz). By the "diagonal sequence trick", there is 

a subsequence, again denoted by p^\ such that pm converges to some pm for 
each m, weakly in £ 3 (M, dz) fl£ p (R, dz) for some 1 < p < 3 and pointwise almost 
everywhere. By Fatou's lemma 



i\z) 



liminf V / dzdz'V m>n (z - z')p m \z)p[ 

i^oo ' * J 
m,n 

> W dzdz'V m , n (z - z')p^\z)p^\z). (3.3) 

m,n 

Observe now that V m < C for some C > and for all m. Moreover, V m (z) < l/\z\, 
so V m G C p for all p > 1. By the weak convergence we can conclude that 

lim / dzV m {z)p m \z) = [ dzV m {z)p^\z) (3.4) 

for each m. Moreover, by the dominated convergence theorem, we get 

HmE / teV m {z)p<$(z) = W dzV m {z)p^\z). (3.5) 

m m 

Using the "diagonal sequence trick" once more, we can use the trace class property 
of the rf's to conclude that there exists a subsequence of and a r DDM = 
(r° DM ) m such that 

T« - TZ DM (3.6) 

in weak operator sense, for each m. It follows from weak convergence that < 
r° DM < 1. Using Fatou's lemma twice, we have 

^ Tr[r D DM] < liminf ^Tr[r«] < N. (3.7) 



in 
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By the same argument 

^Tr[-«9 2 rf DM ] < Uminf J2Tr[-d 2 z r<$]. (3.8) 

m m 

Let now P™ DM denote the density of r° DM . It remains to show that pm = Pm° M 
for each m. From weak convergence it follows that 

(i - #) 1/a rg(i - ^ 2 ) 1/2 - (i - ^ 2 ) 1/2 r™(i - d 2 ) 1 ^ (3.9) 

weakly on the dense set C£°(K). Since the operators are bounded by Tr[(l — 
of)r$] < C, we see that Q holds weakly in £ 2 (M,d;z). With r/ e C£°( M ) con- 
sidered as a multiplication operator, it is easy to see that (1 — 9 2 )" 1 / 2 ?7(1 — 9 2 ) -1 / 2 
is a compact operator (it is even Hilbert-Schmidt). Thus it can be approximated 
in norm by finite rank operators. Using ( |3.9| ) we can therefore conclude that 

.limTr[r^]=Tr[r° DM r ? ], (3.10) 

the sense of distributions. Since we already know that pl2 
converges to p^ pointwise almost everywhere, we conclude that pfif* = p° DM 
for each m. 

We have thus shown that there exists a r DDM with Y2 m Tr[r^ DM ] — N an d 
^,z M [r DDM ] < liminf^^f^frW] = £ DDM . ' ""' □ 

3.2. LEMMA (Uniqueness of the density). The density corresponding to 
the minimizer is unique, i.e., if there are two minimizers and Y^ 2 \ their 
densities pm and pm are equal, for all m. 

Proof. Observe that 

D(p,p) = D(p,p), (3.11) 

where we set p(x) = J2 m l0m( x _i_)| 2 Pm(^)- Using the positive definiteness of the 
Coulomb kernel and the fact that p^ = p^ implies p^m = pm for all m, we 
immediately get the desired result. □ 

Having established the uniqueness of the density, we can now define a lin- 
earized DDM functional by 

^° DM [r] = E Tr ^ DMr -]' ( 3 - 12 ) 

m 

with the one-particle operators 

= - *£ DM (*), (3-13) 

where the potentials are given by 

1%>*(z) = ZV m {z) P™ M {z')V m , n {z - z')dz'. (3.14) 
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3.3. LEMMA (Equivalence of linearized theory). A minimizer r DDM for 
^™ under the constraint ^ m Tr[r m ] < N is also a minimizer for the linearized 
functional (with the same constraint). 

Proof. We proceed essentially as in [|LSY94a|| . For any T 

£»™[T) = £° DM [r] + D( Pr - p DDM , p r - p DDM ) - 5(p DDM , p DDM ). (3.15) 
In particular, for all 5 > 0, 

- £)r DDM + 5T] = (i - 5)£° DM [r DDM ] + ^ n DM [r] - 5( P DDM , P DDM ) 

+5 2 5(p r -p DDM ,p r -p DDM ). (3.16) 

Now if there exists a r<°) with £ TO Tr[r£?] < ^ and ^i? n DM [r (0) ] < ^ n DM [r DDM ] 
we can choose S small enough to conclude that 

££T[(i - ^)r DDM + 5r (0) ] < C M [r DDM ] - 5(p DDM , P DDM ) = ^ M [r DDM ], 

(3.17) 

which contradicts the fact that r DDM minimizes Sg^- D 

We now are able to construct the minimizer of by means of the eigen- 

functions e l m of the one- dimensional operators /i° DM . If p^ < p„ < . . . denote 
the corresponding eigenvalues, there is a p < such that 



F m ~~ ^ ] l e m)( e ml + •^m| e m')( e m'l) (3.18) 



where I m = max{« : p^ < p}, and < A m < 1 is the filling of the last level. 
Since A m is determined by the unique density Pm DM , we immediately get the 
following corollary: 

3.4. COROLLARY (Uniqueness of r DDM ). For any B > and Z > the 

minimizer of the functional ^ M under the condition Yl m Tr[r m ] < N is unique. 

Note that r DDM really depends on the three parameters N, Z and B, but we 
suppress this dependence for the simplicity of the notation. 

Of course, the choice of p is not unique, but I m is unique for every m. In the 
following, we will choose for p the smallest possible, i.e. 

p = max{p^}. (3.19) 

The energy E DDM (N, Z, B) is a convex, non-increasing function in N. Moreover, 
it has the following property. 
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3.5. THEOREM (Differentiability of E DDM ). For N £ N the energy E DDM 
is differentiable in N, and the derivative is given by dE DDM /dN = fi with // given 
above. For N 6 N, the right and left derivatives, dE DDM /dN±, are given by 

Q771DDM a 771DDM 

Proof. For fixed N, Z and B define as above. Let E\\ n (N') be the infimum 

of £jP DM under the constraint E m Tr [ r m] < N ' ■ For < iV_ < iV < N + let Y ± 
be the corresponding minimizers. We have, with < 5 < 1, 

eP™{n + 8(n±-n),z,b) < ^5 M [(i-5)r DDM + 5r ± ] 

= E DDM (N,Z,B) 

+5 (E hn (N±) - E nn (N)) + 0(5 2 ). (3.21) 

Dividing by 5(N± — N) and taking the limit 5^0 followed by N± — > N, we 
conclude that 

dE rm <9E DDM dff DDM < dE lin 

<9iV_ ~ dN_ ' dN + ~ dN+' [ ' ' 

From ( gig ) we infer that E DBM {N',Z,B) > E lin (N') - 5(p DDM , p DDM ) for all 
N', with equality for N = N'. Therefore we get the inverse inequalities in (|3.22| ), 
so equality holds. The assertions of the theorem follow from 

[N] 

E lin (N) = ]T /i(j) + (N- [N]) fi([N] + 1), (3.23) 
i=i 

where [N] is the largest integer < N, and fj,(J) is the j'th element of the set {fi l m } 
in increasing order. □ 

Corresponding to the DDM minimizer we define the three-dimensional DDM 
density as 

p DDM (x) = ^|^( Xi )|Vr(4 (3-24) 

m>0 

3.6. THEOREM (Ordering of the Assume /i° DM has at least M eigen- 
values. Then, for m > 1, 

M _ M r 

+ (™ + - (2^ + 1)a4) < - ^ E / ^ DM |^<J 2 ( 3 - 25 ) 

8=1 i=l ^ 

i/ie understanding that \i % m _ x (f4n+i) = ® ^m~t (^m+i) ^ as ^ ess th an i 
eigenvalues). The analogous inequality for m = is 

m m _ m r 

E (/4 - /4) < Z E l e o(°)! 2 - f E / P DDM |^eo| 2 . (3-26) 

i=l m=l i=l 
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Proof. With A(2) the two-dimensional Laplacian, one easily computes that 

(2 J B)- 1 A (2) |0 m | 2 = m|0 m _x| 2 + (m + l)|0 m +i| 2 - (2m+ l)|0 m | 2 (3.27) 

for m > 1. Multiplying ( |3.27| ) with $(x) = Z/|x| — p DDM * |x| _1 and integrating 
over xj_, we therefore get, for any density matrix 7 (recall the definitions (|1.14j) , 

qnsp, <$n$ and o), 



Tr [(mh™? + (m + - (2m + 1)/£ DM ) 7] 



-(2B)~ l J dx$(x)p 7 (^)A (2) |0 m (x ± )| 2 

-| / dxp DDM |0 m | 2 p 7 + ^ / d ^° M J^, (3-28) 



where we used partial integration for the second step, and the fact that <f) m (0) = 
for m > 1. To treat the last term in ( 3.28|) , note that the function 



M 

w — > 

i=l 



M „ 

£ y d^ DM W|e^(^ + ^)r (3.29) 



has its maximum at w — 0, because otherwise one could lower the energy by 



shifting the eigenvectors e l m . Therefore the second derivative of (|3.29|) at w = 



is negative. Setting 7 = J2 i=1 \ e m)( e m\ we see that the l as t term in (|3.28| ) is 
negative, so we can conclude that 

Tr [{mh^t + (m+ 1)/^ - (2m + 1)/£ DM ) 7] < ~ J p DDM |0 m | 2 p 7 - 

(3.30) 

By the variational principle (|3.25|) holds. The proof of (|3.26|) is analogous, con- 



sidering also the contribution from 0o(O) in Q3.28|) . □ 

As a corollary, we immediately get 
3.7. COROLLARY (Monotonicity of YaU /4 in m )- YaU 

fj, l m is increas- 
ing in m for all M, and strictly increasing as long as it is < 0. Moreover, 
t ddm = Q f orm >N. 



Proof. This follows immediately from Theorem 1370, noting that linim^oo [i l m = 



for all z. □ 

3.8. Remark (Concavity of Y^=itAn)' The result ( |3.25| ) (together with the 
monotonicity) implies concavity of y} m in m. More precisely, 

m m + 1 1 1 1 1 1 / \ 

2mTl /im - 1 + 2mTl /im+1 + 2(2m + 1) (/Xm " 1 " ^ m+l) ~ ^ (3 ' 31) 
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(and the same holds for ^2 i=1 ^ m )- This is the analogue of Prop. 2.3 in ||LSY94a 
which states that — /ii(xj_) is a increasing and concave function of |xj_| (note the 
different sign convention). 

We now introduce the parameters A = N/Z and rj = B/Z 3 . The next theorem 
deals with the 77 — > 00 limit of E DDM with fixed A. To prove it we need the 
following Lemma. 

3.9. LEMMA (Convergence to delta function). Let L = L(rj) be the solu- 
tion of the equation 

rfl 2 = L(rj) sinh(L(r7)/2). (3.32) 

Let ip G H^^dz), with X = J \ip\ 2 and T = J \dip/dz\ 2 . Then, for all m > 0, 
Z > andr] = B/Z 3 , 

l^(0)| 2 -^ / V m (z/LZ)\iP(z)\ 2 dz 



Proof. After an appropriate scaling, this is a direct consequence of PSYOO 



Lemma 2.1, using the estimates 



IxJIVxI-dx, < f-^^ (3.34) 



and 

1)1/4 



J \M^)\ 2 \*±\ 1/2 dx± < 2 {m ^} . (3.35) 



□ 



3.10. THEOREM (The limit 73 > Z 3 ). For all A > 

lM (AZ, Z, 
r?^oo Z 3 (In 77) 2 

uniformly in Z . 



lim ^(AZ.Z,*) = (3M) 



Remark: The uniformity in Z will be important for the proof of Theorem p. 12 
It is non-trivial in contrast to DM, where one has the scaling relation ( |2.8|) , 
which implies that the left hand side of ( |3.36| ) (with DDM replaced by DM) is 
independent of Z. 
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Proof. The lower bound is quite easy, using the results of ||LSY94a |. As shown 
in Section El we have 



jj/ DDM (AZ, Z, r/Z 3 ) > E™(\Z,Z, V Z 3 ) £DM ^ 



Z 3 



z 3 



(3.37) 



where we have used the scaling properties of E DM . It is shown in ||LSY94a l that 
the right hand side of ( |3.37| ) divided by (liirj) 2 converges to E ns (X). 

For the upper bound we assume N £ N for the moment. We use as trial 
density matrices 



r m 0, z') 
r m (z, z') 



L^p HS (LZz)^p ns (LZz'), < m < N - 1 



N 

0, m > N 



(3.38) 



where p is given in ( |2.10|) and L = L(rj) is defined in ( |3.32| ). The kinetic energy 
is easily computed to be 

Z 3 L 2 j \d^/p^/dz\ 2 . (3.39) 
For the attraction term we use Lemma |3.9| to estimate 

> p hs (o)-§ 

for some constant C' x depending on A. 

For the repulsion term we first estimate 



(m + 1) 



1/4 



(3.40) 



/5\ 2 y e(^2N/B - |x ± |)e( v ^iV7B - |x'J) 



dx i dx' 



±5 



0<m,n<N-l 

(3-41) 

which follows from monotonicity of l/|x| in |x_|_| and the fact that Y2 m l^ml 2 < 
B/2n. Therefore we have 



Z^j2 j V m)n (z - z')p m (z)p m (z')dzdz' 

0<m,n<N-l 

fU(z-z'))p HS (z)p us (z')dzdz', 
where we set £ = (77/ A) 1//2 j L, and the function / is given by 



(3.42) 



f(z) = (27T) 



-2 



e(V2- 


x± |)e(V2- 


x'J) 


\/l x ^ - x ± 


2 + z 2 



dx_|_dx'|. 



(3.43) 
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We now claim that L 1 £f(£,z) — > 5(z) as 77 — > oo. Since f(z) < l/\z\ we have, 
for any x £ H 1 H 

^ y f(^z)x(z)dz- X (0)L- 1 J ^f(z)dz 

= L- l i f f^z)( X (z) - X (0))dz + L"^ / /(£z)x(*)d* 
</-i "'|z|>i 

< const. L -1 ^ / f(£z)z x/2 &z + L- x f \x{z)\dz < const./,- 1 . (3.44) 

7-1 ■/|2|>l 

Now Ji z |<^/ 1 as 77 — > oo, which proves our claim. And since £(77) ~ In 77 
for large 77, this finishes the proof of Theorem [3.10| , in the case where N is an 
integer. The proof for N ^ N is analogous, using ( |3.38|) with the density matrix 
corresponding to m — [N] multiplied by iV — [N] as trial density matrices. □ 



3.11. COROLLARY (HS limit of the density). For fixed A = N/Z 



ri'-^oo Z 2 In 77 



In n 



,DDM, 



z/Zlni]) = p 



HS/ 



(3.45) 



in the weak C l sense, uniformly in Z . 



Proof. The convergence of the densities in ( |3.45| ) follows from the convergence 
of the energies in a standard way by considering perturbations of the external 
potential (cf. e.g. [|LSY94aj| ). Moreover, since the convergence in ( |3.36| ) is uniform 
in Z, ( |3.45| ) holds for any function Z = Z{0), so we can conclude that (ft.45|) holds 
uniformly in Z, too. □ 

Using the results above we can now prove the analogue of Theorem 4.6 in 
SY94 £ 



3.12. THEOREM (r° DM has rank at most 1 for large 77). There exists a 
constant C such that 77 > C implies that T^ DM has rank at most 1 for all m. 



Proof. We first treat the case A < 2. From Theorems |3.5| and |3.10| and the fact 
that E DDM (XZ, Z,rjZ 3 ) is convex in A, we get 



// 



i^o Z 2 (In 77) 2 



dE Uh (\) 
dX 



< 0. 



(3.46) 



Suppose that \i is not the ground state energy of some /i° DM - Then \i > —Z 2 /4, 
because the second lowest eigenvalue of h l ^P M is equal to the ground state energy 
of the three-dimensional operator 



- A-$ 



DDM 

m 



(W). 



(3.47) 
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This follows because $^ (z) is reflexion symmetric, so the eigenvector corre- 
sponding to the second lowest eigenvalue, u m (z), has a node at z = 0. Therefore 
u m (|x|)/|x| is an eigenvector of ( |3.47|) , and because it does not change sign, it 



must be a ground state. Since $° DM (|x|) < Z/|x|, the ground state energy of 
( |3.47| ) is greater than — Z 2 /4. So /i/Z 2 (ln rf) 2 would go to zero as r\ — > oo, in 



contradiction to ( |3.46| ). Therefore there exists a constant C such that rj > C 



implies the assertion to the theorem. This constant can be chosen independent 
of Z, because the limit ( |3.46| ) is uniform in Z (by the same argument as in the 
proof of Corollary |3.11|) . 

Now assume that A = 1 + 5 for some 5 > 0. From Corollary |3.11| we infer that 
for large enough 77 there is some cs > such that 

E/ P° DM >(1 + ^)^ (3-48) 

where 5 = min{l, 5}. We now will show that for rj large enough /i„ DM has at most 
one eigenvalue, for all < m < N. By the same argument as above, we need 
to show that the three-dimensional operator ( |3.47| ) has no eigenvalues. Using 
Kn( ;Z ) — 1/1 2 1 an d (cf- the next section) 

V m Jz) > 4=Kn+n(^/v / 2) > — = = 1 ===== > 



V2 ^/z 2 + 2{m + n + l)B- 1 ~ y^z 2 + 2(2N + ^B' 1 

(3.49) 

we can use fl3.48|) to estimate 



JJL».\i,i^i, / Z Z(l+ gtf) 



< UM (|x|) < ^ - , V 2 ' (3.50) 

11 ; (|x| + csiZlnri)' 1 ) 2 + 2(2N + l)^- 1 

Therefore, for 77 large enough, 

*r*(i*i)<{ z/ J x| Ul^ (Zh, " ) " 1 (3.5D 



By means of the Cwikel-Lieb-Rosenbljum bound [RS7S] we can estimate the 
number of negative eigenvalues of (|3.47| ) as 

const. J |$° DM (|x|)|fdx < c' s (\n V y 3 / 2 , (3.52) 

which is less than 1 for 77 large enough. □ 



3.13. Remark (Chemical potential for large B/Z 3 ). The theorem above, 
together with Corollary |3.7| , shows that for B/Z 3 large enough, the chemical 
potential is given by the ground state energy of /i<° r ^i > , i.e. \i = /i<jv-i>) where 
< N > denotes the smallest integer > N. 
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4 Maximal negative ionization 

The DDM energy is convex and monotonously decreasing in N. Because D(p,p) 
is strictly convex in p, _£; DDM is strictly convex up to some iV = N C (Z,B), and 
constant for iV > iV c . By uniqueness of r DDM the minimizer for N > N c is equal 
to the one with N = N c . In particular, ]T m Tr[r° DM ] = min{iV, iV c }. 

The "critical" N c measures the maximal particle number that can be bound 
to the nucleus. We will proceed essentially as in ||S00|| and use Lieb's strategy 



84 1 to get an upper bound on N c . In addition, the following Lemma is needed. 
Throughout, we use various properties of V m stated in ||BRW99|| , Sect. 4., and 
proven in ||RWO0|| . 



4.1. LEMMA (Comparison of V m and V m , n )- 



1 + 77T^ + 77^ + 77T77 ) V m Jz-z')>\ (4.1) 



V m {z) V n {z) V rn (z') V n {z 
Proof. Using the definition of V myn it can be shown (cf. ||P82|| ) that 



ra+n 



V2V mtn {V2z) = C M Z ) ( 4 - 2 ) 



i=0 



for some coefficients c« > that fulfill ^ q = 1. Since V m > V m+ i for all m 
( ||BRW99|| , 4b) we get ^ 

V m , n {z) > ^=V m+n {z/V2) > W m+n {z/2), (4.3) 

where we used the fact that aV m (az) < V m (z) if a < 1 ( PRW99]] , 4g). Moreover, 
using convexity of l/V m+n ( ||BRW99| , 4i), we arrive at 

V m , n (z - Z') > {V m+n {z)~ l + Vrn+niz')- 1 )' 1 ■ (4.4) 



The assertion ([4.1| ) follows if we can show that 



< + (4.5) 



V m+n (z) V m (z) V n (z) 

This is of course trivial if n or m equals zero. If n, m > 1 we use V z 2 + m > 
V m (z) > Vz 2 + m + l^ ( PRW99|| , 4a) to estimate 



^m+n (^0 



< Vz 2 + m + n + l<2Jz 2 + (y/m + y/n) 2 /A 



1 1 



< Vz 2 + m + Vz 2 + n < — — + — — , (4.6) 

V m (z) V n (z) 

which finishes the proof. □ 
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4.2. THEOREM (Critical particle number). 

1 dE DDM (N c ,Z,B) 



Z < N r < AZ — 



Nr. 



dZ 



(4.7) 



Remark: The factor 4 stems from the symmetrization of ( 4.1|) in m and n. Due 



to this symmetrization one could expect that Lemma [4.1| holds with 1 replaced 
by 2 on the right hand side. This would imply that 4Z could be replaced by 2Z 
in 



(4.t 



Proof. Let ej„ denote the eigenvectors of /i° DM , i.e. 



^DDM . = 



Multiplying ( |4.8|) with e^/V^ and integrating we get 
1 



2>E 



1 



4(W - z'K DM (z')dzdz' + (ej — (-^)|0, (4.9) 



-it- / \ 77t \ / • - irn \ I 1 \~mi-rr 

*rn\Z) *m 



where we used that fi l m < 0. Since 1/V m is convex and |z|V^(z) — > 1 as |z| — > oo 
we have |(l/V^i)'| < 1. Using this and partial integration we can estimate the 
last term in (|4.9|) by 



V, 



VL 



Vrr. 



I m m I 



> -- i v m ( z yjzYdz. 



(4.10) 



Summing over all m and i (and, according to Eq. ( |3.18| ), multiplying the factor 
corresponding to the largest p l m by A m ), we arrive at 



m,n m m 

(4.11) 

Note that the last term in ( 4.11 ) is equal to dE DDM /dZ. To treat the first term 
in ( 4.11 ) we use symmetry and Lemma |4.1| to get 



E 



i\ „DDM/ J\ 



V m z) 



1111 
+ tt~i — r + „ , + 



V m (z) V n (z) V m (z') V n (z>) 



xp» DM (z)V m>n (z - z')p» DM (z')dzdz' > - A N 2 . (4.12) 
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Inserting this into ( 4.11|) and dividing by N/A we arrive at 



1 <9K DDM 

N ^ 4Z ~ N— ■ < 4 - 13 > 

The lower bound on N c is quite easy. We just have to show that ^j^ffi has a 
bound state if N < Z. Using if)(z) = exp(— a\z\) with a > as a trial vector we 
compute 

= a-Z f V [N+1] (z)e- 2a \ z \dz + J2 [ V n , [N+1] (z - z')p» BM (z')e-^dzdz' 
<a-Z J V [N+1] (z)e- 2alzl dz + Nmax J V ni[N+1] (z)e- 2a ^dz. (4.14) 

Since 

lim — 1- / V [N+1] e~ 2a ^ = lim — 1- / V nAN+1] e~ 2a ^ = 1, (4.15) 
o-»oln(l/a)j a->oln(l/a) 7 

(VH^w+iil^) wm be negative for small enough a, if iV < Z. □ 

4.3. Remark (Explicit bound on N c ). The concavity of E BBM in Z implies 
that dE DDM (N,Z,B)/dZ > E DDM (N, 2Z, B)/Z. Using ( gig) and the bounds 
on £' DM given in |[LSY94a(| , Thm. 4.8, (|4.7p implies the upper bound 

iV c < AZ (1 + Cmin {(B/Z 3 ) 2/5 , 1 + [\n(B /Z 3 )] 2 }) (4.16) 

for some constant C independent of Z and B. 



4.4. Remark (Upper bound on N BU ). The upper bound ( 4.16|) holds also 



for N® M , the critical particle number in the DM theory. In fact, the convergence 
in Q2.13| ) implies that 

Nr<Z L m M^n) (4.17) 

for all fixed rj = B/Z 3 . 

5 Upper bound to the QM energy 

We now show that E DDM is an upper bound to the quantum mechanical ground 
state energy. In fact it is even an upper bound to -E^nf By Lieb's variational 
principle |[L81 | 



Ef> < Tr[(H A - Z|x|- 1 ) 7 ] + \ J 7(x ; x) 7 (x^)-| 7 (x,xQP dxdx/ ^ 
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for all density matrices < 7 < 1 with Trpy] < N. We choose 

7(x, x') = ^ 0m(x_ L )0 m (x'_ L )r m (z, z'). (5.2) 

m 

Since n 7n = 7, we get an upper bound even for £^ nf . Omitting the negative 
"exchange term", we compute 

^c Q onf(^^5)<^ M [ r ]' (5.3) 

where T = (T m ) m . Note that Tr[j] = 5] m Tr[r ro ]. Therefore we immediately 
conclude that 

£ c Q onf (N, Z, B) < E DDM (N, Z, B). (5.4) 

6 Lower bound to the QM energy 

To get a lower bound on the QM energy we need to estimate the two-body 
interaction potential in terms of one-body potentials. One way to do this is to 
use the Lieb-Oxford inequality |LQ81| together with the positive defmiteness of 
the Coulomb kernel: 

<*| ^Ix.-x.r 1 ^) >2D(p,pv)-D(p,p)-1.68 I pf, (6.1) 

i<j J 

where we can choose p = p DDM . However, if \1/ is an (approximate) ground state 
wave function in the lowest Landau band, the error term J p^ 3 will in general 
be greater than the energy itself. More precisely, if 7^ denotes the one-particle 
density matrix of we can estimate 

*» < (/^f(/-f 

/ o \ 1/6 

< KS 2 Tr[-^]J iV 5 / 6 

< const.Z 1 / 5 ^ 14 / 15 ^ 2 / 5 , (6.2) 



where we used |[LSY94a|| , Lemma 4.2, to estimate j p| and the kinetic energy. 



For large enough B this bound is of no use. 

That is why we follow the method of |[LS Y9 lajj to get an improved bound for 



large B. Their result is that for wave functions \1/ that satisfy (^f\H^f) < the 
bound 

(*| I*-**!" 1 *) > 2^(P DM ,P*) - D(p™, P ™) -c x Z*' a (1 + (In??) 2 ) (6.3) 

i<j 
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holds. The proof of this result uses only properties of p DM which hold also for 
p DDM . Therefore ( |6.3|) holds also with p DM replaced by p DDM ; possibly with a 
different constant. So we can estimate 

E? oni (N,Z,B) > inf^l^^-Z/W+p^Hxr 1 )*) 

-D(p»™,p»™)-R L , (6.4) 
where the infimum is over all \1/ with Uq^/ = \1/ and \\^\\ = 1, and 

R L = c x min{Z 17 / 15 5 2 / 5 , Z 8 / 3 (l + (In r/) 2 ) }. (6.5) 

Because — Z/\x\ + p DDM * |x| -1 is a one-particle operator that is invariant 
under rotations around the z-axis, we can restrict ourselves to considering Slater 
determinants of angular momentum eigenfunctions, which leads to 



E? oni (N,Z,B) > MS^ M [T]-D(p^ M ,p^ M )-R L 

iDDM / 



= E UUM (N,Z,B) -R L . (6.6) 

6.1. Remark (Magnitude of the "exchange term"). One may ask for the 

optimal magnitude of Rl such that a bound of the form ( |6.6|) is valid. Rl is an 
upper bound on the difference between E® oni and E DDM , which is given by the 
exchange energy 

h = D{py,pv) - (tf | l x * ~ Xj| _1 *>. ( 6 - 7 ) 

i<j 

where \I/ is an (approximate) ground state wave function of H. The exchange 
energy is roughly N times the self energy of the charge distribution of one particle. 
This distribution has the shape of a cylinder with diameter R ~ B^ 1 ^ 2 and length 
L ~ max{Z- 2 / 5 J B' 1 / 5 ,Z- 1 ln( J B/Z 3 )- 1 }, if we assume N ~ Z (see also ||LSY95 



for heuristic arguments). Note that L ^> R is equivalent to B ^> Z 4//3 . In this 
case, we expect an exchange energy of the order 

ZL-HniL/R) ~ Z 7/5 B 1/5 \n{B/Z 4/3 ) for Z 4/3 < 5 < Z 3 

ZZT 1 ln(L/.R) ~ ^ln^/Z^ln^/Z 2 ) for 5 > Z 3 . (6.8) 

Hence the exchange energy should be of order Z 1 ^ 2 (E DDM ) 1 ^ 2 x some factor loga- 
rithmic in B. We conjecture that (at least for appropriate \1/ in the lowest Landau 
band) 

Iq, < const.— pq,x some factor logarithmic in B, (6-9) 



5 , 

which is precisely of the correct order ( |6.8| ). This is in accordance with results 
on the homogeneous electron gas in a magnetic field | pG71| , |FGPY92 |. 
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